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Abstract. According to Einstein, the notions of geodesic, parallel transport (afhnc 
connection), and curvature of space-time manifold have a pure geometric origin and 
do not correlate with any electromagnetic concepts. At the same time, curvature is 
generated by matter which is not affiliated with the spacetime geometric concepts. For 
this reason, the fundamental constant c entering the geometric and matter sectors of 
general theory of relativity have different conceptual meanings. Specifically, the letter 
c in the left side of the Einstein equations (geometric sector) entering the Christoffcl 
symbols and its time derivatives is the ultimate speed of gravity characterizing the 
upper limit on the speed of its propagation as well as the maximal rate of change of 
time derivatives of the metric tensor, that is gravitational field. The letter c in the right 
side of the Einstein equations (matter sector) is the maximal speed of propagation of 
any other field rather than gravity. Einstein's general principle of relativity extends 
his principle of special relativity and equates numerical value of the ultimate speed 
of gravity to that of the speed of light in special theory of relativity but this general 
principle must be tested experimentally. To this end we work out the speed of gravity 
parameterization of the Einstein equations (c g -parameterization) to keep track of the 
time-dependent effects associated with the geometric sector of general relativity and 
to separate them from the time-dependent effects of the matter sector. Parameterized 
post-Newtonian (PPN) approximation of the Einstein equations is derived in order to 
explain the gravitational physics of the jovian deflection VLBI experiment conducted 
on September 8, 2002. The post-Newtonian series expansion in the c g -parameterized 
general relativity is with respect to a small parameter that is proportional to the ratio of 
the characteristic velocity of the bodies to the speed of propagation of the gravitational 
interaction c g . The Einstein equations are solved in terms of the Lienard-Wiechert 
tensor potentials which are used for integrating the light-ray propagation equations. 
An exact analytic expression for the relativistic time delay in the propagation of a 
radio wave from the quasar to an observer is calculated under the assumption that the 
light-ray deflecting bodies move with constant velocities. A post-Newtonian expansion 
of the time delay proves that in general relativity the time delay is affected by the speed 
of gravity already to the first order in l/c g beyond the leading (static) Shapiro term. 
We conclude that recent measurements of the propagation of the quasar's radio signal 
past Jupiter are directly sensitive to the time-dependent effect from the geometric 
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sector of general relativity which is proportional to the speed of propagation of gravity 
c g but not the speed of light. It provides a first confirmative measurement of the 
fundamental speed c of the Einstein general principle of relativity for gravitational 
field. A comparative analysis of our formulation with the alternative interpretations 
of the experiment given by other authors is provided. 
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1. Introduction 



The relativistic VLBI experiment for ultra-precise measurement of the deflection of light 
caused by the time-dependent gravitational field of moving Jupiter was conducted on 
September 8, 2002 by the National Radio Astronomical Observatory (USA) and the 
Max Plank Institute for Radio Astronomy (Germany) pp. The idea of the experiment 
was proposed by Kopeikin [2] who noted that a moving gravitating body interacts 
with a light particle (photon) not instanteneously but with retardation due to the 
finite speed of propagation of the body's gravitational field to the light particle (see 
Fig. fJJ and Fig. |2J). Thus, two null (causality) cones are essential in the problem 
of the calculation of relativisitc deflection of light by moving bodies - the null cone 
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associated with propagation of light (light null cone) and the null cone associated with 
propagation of gravity (gravity null cone) from the light-ray deflecting body (see Fig. 
EJ). The retarded position of the light-ray deflecting body from which it deflects light, is 
connected to the present position of the light particle via an equation of the gravity null 
cone which is defined by the causal solution of the wave Einstein's equations taken in 
the form of the Lienard-Wiechert potentials. The retarded position of the body is taken 
on its orbit at the retarded time associated via Lienard-Wiechert retarded time equation 
with the finite speed of gravity c g which must be numerically equal in general relativity 
to the speed of light c in accordance to the Einstein general principle of relativity ||. 
The goal of the experiment was to confirm the Einstein general principle of relativity 
by making use of the close celestial alignment of Jupiter and the quasar J0842+1835 
and to prove that in accordance to this principle the ultimate speed of propagation of 
gravity is the same as the speed of light. This paper discusses the gravitational physics 
of the experiment in more detail by calculating exact expression for the VLBI time 
delay in terms of the speed of gravity parameter c g that is measured with respect to the 
speed of light c which numerical value is well-known from laboratory experiments. The 
post-Newtonian expansion of the VLBI time delay explicitly in terms of the expansion 
parameter e = c/c g leads to a more profound understanding of our experimental results 

m 

The paper is organized as follows. In section 2 we discuss the speed-of-gravity 
parameterization of the Einstein equations that is used in this section for derivation of 
the linearized gravitational field equations, the metric tensor, and the affine connection. 
We formulate the law of the parallel transport for test particles and derive the exact 
formula for the gravitational time delay of light in section 3 in terms of the present 
and retarded variables. In section 4 we focus on its post-Newtonian expansion and 
dependence on the speed-of-gravity parameter c g . Then, the differential VLBI time delay 
between two VLBI stations is derived and we prove that our original Lorentz-invariant 
formulation of the VLBI delay ,2] is identical with the post-Newtonian presentation of 
this paper. General relativistic interpretation of the experiment associated with the 
speed of gravity is discussed in section 5. Non-general relativistic interpretation of the 
experiment in two-parametric model of gravity is outlined in section 6. Alternative 
interpretations of the experiment given by other authors are analyzed in section 7. 

1.1. Fundamental Speed c and Its Classification 

A correct and unambigious interpretation of the jovian deflection experiment is 
inconceivable without peer study of the nature of the fundamental constant c appearing 
in various differential equations of theoretical physics. Having one and the same name 
this constant has different physical meanings and, in fact, may tacitly represent a set of 

|| In general relativity, both the speed of gravity and the speed of light are denoted by the same letter 
c. To keep track of the effects associated with the speed of gravity and to avoid confusion we shall 
denote the speed of gravity by the symbol c g as contrasted to the symbol of the speed of light c. 
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fundamental constants of nature characterizing various physical properties of the world 
in which we live. It is vitally important to carefully distinguish different facets of c from 
each other. 

Recently, Ellis and Uzan j^j have suggested a classification scheme in order 
to prevent confusion between different c that would facilitate a correct physical 
interpretation of experiments in fundamental physics. Extending the idea of Ellis and 
Uzan [3J we define a set of c in an electromagnetic field theory 

(i) c ew - the speed of electromagnetic waves in vacuum , 

(ii) q - the ultimate speed of light (electrodynamic constant) , 

(iii) cm - the coupling constant between electromagnetic field and electric current , 

and a set of c in a gravitational field theory 

(i) - the speed of gravitational waves in vacuum 

(ii) c g - the ultimate speed of gravity (gravitodynamic constant) , 

(iii) ce - the coupling constant between gravitational field and matter. 

Physical meaning of the constants q and c g is that they determine the maximal rate of 
change of electromagnetic and gravitational fields respectively. It means that q appears 
in front of the time derivatives of electric and/or magnetic fields while c g will appear in 
front of the time derivatives of the gravitational field (the metric tensor). 

The most general theory of electromagnetism might, in principle, have three 
constants: c ew , q, and Cm, to be different. Maxwell's theory of electromagnetism 
establishes an exact relationship: c ew = q = Qvi, which has been confirmed with an 
unparalled degree of precision in a multitude of laboratory experiments. For this reason, 
we accept that Maxwell's theory has a single constant, which name is the speed of light 
c. We emphasize that the speed of light characterizes all electromagnetic phenomena 
but not only the propagation of electromagnetic waves in vacuum. 

Before Maxwell created his theory the nature of light was not known but its speed 
was measured fairly well in 1676 by a Danish astronomer, Ole Romer, working at the 
Paris Observatory It is clear now that what Romer had measured was c ew . This was 
first indicated by Maxwell who noted that the product of the electric permittivity e and 
magnetic permeability no of vacuum together determine the velocity which is in excellent 
numerical agreement with the speed of light. This led Maxwell to the conclusion that 
light is electromagnetic wave. According to our notation the product (e fi Q = q and 
it might be different from c ew in a different world. In our world, however, c ew = Q = c. 

It is likely that the most general theory of gravity should, in principle, also contain 
three constants: c gw , c g , and ce to be different. Einstein's theory of general relativity 
assumes that c gw = c g = ce and equates c g and c. These assumptions do not contradict 
timing observations of binary pulsars with the accuracy approaching 0.4% jjjj. However, 
it does not mean that pulsar timing measures the numerical value for Cgw with this 
precision because the timing model of binary pulsar observations do not incorporate the 
constants Cgw, c g , and ce explicitly. Much more work is required to incorporate these 
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parameters in the processing software of timing observations in order to measure them 
separately. 

We shall continue our analysis of the VLBI experiment under discussion in the 
framework of general relativity by parameterizing its equations with a single parameter 
c g = c gw = c g = ce- We notice that the ratio G/c 2 in the Einstein equation is fixed 
by the principle of correspondence of the Einstein equations with Newtonian gravity. 
Therefore, the parameter c g can appear either in the left side of the Einstein equations 
as a coefficient in time derivatives of the metric tensor or in the right side of the Einstein 
equations in terms proportional to G/(c 2 c g ) and/or G/(c 2 c 2 ) (see [2] and next sections 
for more detail). 

We call the parameter c g as the (ultimate) speed of gravity in close analogy with 
the electromagnetic theory, where a similar unifying constant c is called the speed 
of light. We emphasize that the speed of gravity c g pertains to all time-dependent 
gravitational phenomena but not only to the propagation of gravitational waves in 
vacuum. Any theory of gravitational field in which c gw 7^ c g 7^ Ce will have equations 
with mathematical properties different from general relativity. At any theory Einstein's 
principle of relativity demands c gw < c g so that those theories for which this inequality 
is violated must be considered as invalid. We shall briefly discuss interpretation of the 
jovian deflection experiment in two-parametric model of gravity in section 13 



1.2. The Einstein Equations and the Speed of Gravity Concept 
Einstein equations have the symbolic form 

GWc] = ^-T, v [c) , (1.1) 

where G^ v is the Einstein tensor, T^ u is the stress-energy tensor of matter, G is the 
universal gravitational constant, and c is a fundamental speed that is numerically equal 
to the speed of light in vacuum. The Einstein tensor is a (non-linear) differential operator 
acting on the metric tensor g^ u and its first and second derivatives with respect to time 
and spatial coordinates. It represents a geometric sector of the Einstein equations and 
can be written symbolically as 

{A d 2 B d 2 d 2 1 

#W + iow^ + c -d^f g ^ x) + non - linear terms ' (L2) 

where A, B iy are constants depending on the choice of gauge conditions imposed 
on the metric tensor. Einstein's theory of general relativity presumes that all time 
derivatives 8q of the metric tensor in the Einstein equations are coupled with c, that is 
do = c~ 1 dt as shown in equation p. 21 The ongoing convention is to call the constant 
c "the speed of light" 3j. However, the metric tensor g^ v is not simply a geometrical 
object but represents one of the most fundamental objects in physics - the gravitational 
field. For this reason, the constant c in the Einstein tensor characterises the speed of 
gravitational field and has to be associated with the speed of gravity rather than with 
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the speed of light which has an electromagnetic nature and is physically irrelevant for 
the Einstein tensor. 

On the other hand, the stress-energy tensor T^ v is defined locally as a special 
relativistic object and can not physically depend on the speed of gravity in a direct way 
because gravitational field is not localized. Nonetheless, T^ u can depend on the speed of 
gravity indirectly through the metric tensor g^ v . This dependence may be important in 
higher orders of the post-Newtonian approximation scheme. Thus, we have to keep in 
mind that the fundamental speed c entering the (special relativistic) definition of T^ u is 
the speed of light. Confusion in the interpretation of physical effects can arise, however, 
if one keeps the same notation for the speed of gravity and the speed of light. In order 
to avoid it, one will denote the speed of gravity as c g and the speed of light as c in 
accordance with convention adopted in section fTTTI Then, the Einstein equations (jl.lj) 
assume the symbolic form 

where we have explicitly shown the presence of the speed of gravity c g in the Einstein 
tensor, which must be used in the time derivatives of the metric tensor, that is 
do — > So = c~ 1 d t , so that equation (jl.2j) is recast in the following form 

{A d 2 B d 2 d 2 1 

^ + + ^g-g- j M*> *) + ^n-Hne^r terms . (1.4) 

Of course, the geometrical formulation of the general theory of relativity and the Einstein 
equations are still valid if one keeps the numerical values for c g and c the same. 

A problem will arise in developing c g -parameterization of general relativity with 
the speed of gravity c g taken as a parameter running from c g = oo (Newton) to c g = c 
(Einstein). Sometimes, the Newtonian limit of general relativity is erroneously stated 
as if the speed of light c —>■ oo but this is incorrect. One has to recognize that the 
Newtonian limit of general relativity is, in fact, the limit when c g — > oo while the speed 
of light c is held constant. Indeed, gravity is a physical phenomena characterizing global 
properties of the space-time manifold while the locally flat Minkowski geometry exists at 
each point of the manifold irrespectively of whether gravity propagates instanteneously 
or with finite speed. A rather simple argument supporting this observation can be taken 
from classical astronomy which is based on the Newtonian theory of gravity (c g = oo) 
that does not contradict to existence of such phenomena as the aberration of light - the 
effect confirming that light has a finite speed. 

Development of the c g -parameterization of general relativity demands that the post- 
Newtonian expansion of the Einstein equations with the parameter e = c/c g running 
from e = (Newton) to e = 1 (Einstein) must preserve all the geometric properties of 
general relativity - the Bianchi identity, the relationship between the Christoffel symbols 
and the metric tensor, the gauge and coordinate invariance of the left side of the gravity 
field equations. To preserve the law of conservation of the stress-energy tensor one must 
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introduce the c g parameter to the right side of the Einstein equations explicitly jB] 

G A c z\ = ^-©AB/tCg] , (1-5) 

where ©^[cg] is the stress-energy tensor of matter parameterized by c g in order to make 
the law of conservation for Ml/ consistent with the Bianchi identity for the Einstein 
tensor G M ^[c g ]. This parameterization scheme is explained in the next sections in more 
details. 

Our approach is an extention of the fiber-bundle parameterized approach to general 
relativity developed earlier by other authors (see, for instance, [HIE])- I n this approach a 
parameterized sequence of solutions of Einstein's equations is vizualized as a fiber bundle 
whose base space is the real line (parameter e) and whose fibers are diffeomorphic to 
R 4 , each being space-time for a particular value of e. The fiber e = is Minkowski space 
with a (non-degenerated) Newtonian limit. The fiber e = 1 is a space-time manifold 
of general relativity we are interested in. It is simple to identify physical origin of 
relativistic effects connected with the finite speed of propagation of gravity by tracing 
terms marked by the parameter e. 

1.3. Statement of the Problem and the Lienard-Wiechert Solution of the Einstein 
Equations 

The experimental problem of measuring the relativistic deflection of light in the jovian 
deflection experiment is formulated as follows (see FigHJ). Light rays are emitted by a 
quasar (QSO J0842+1835) at the time t and move to the network of VLBI stations 
located on the Earth. As the light moves it passes through the variable gravitational 
field of the solar system (Jupiter, Sun, etc.) and is received by the first and second VLBI 
stations at the times t\ and ti respectively. The gravitational field of the solar system 
causes a delay in the propagation of radio signals - the effect discovered by Shapiro [Oj. 
We noted that the present-day accuracy of phase-reference VLBI measurements is good 
enough to detect a relativistic correction to the Shapiro time delay depending on the 
retarded orbital position of Jupiter caused by the finite speed of propagation of gravity 
We calculated this correction [2j for Jupiter by making use of the retarded Lienard- 
Wiechert type solution of the Einstein gravity field equations (see equation ()2.23|) in 
section l2~3*j) for the perturbation h^ v of the metric tensor. This perturbation caused by 
a massive point-like body a at the field point (t, x) is 

x ) ~ — -i , x 7T > (1-6) 

c 4 r a - c g v a (s) ■ r a (s) 

where Tfi) is a stress-energy tensor taken on the world line of the body, r a (s) = x—x a (s), 
r a (s) = \r a (s)\, x a (s) and v a (s) are coordinates and velocity of the body taken at the 
retarded time 

s = t- —\x - x a (s)\ . (1.7) 

Crr 
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The retardation in this equation is caused by the finite value of the speed of propagation 
of gravity, c g , in accordance with the causal nature of the Lienard-Wiechert solution of 
the Einstein equations. 

In the general theory of relativity the speed of gravity c g is postulated numerically 
equal to the speed of light c in vacuum, i.e. c g = c. For this reason, Einstein 
used everywhere in his equations only the symbol c. Experimental gravitational 
physics, however, uses light or radio waves for measuring various characteristics of the 
gravitational field. The Einstein convention may be misleading if one wants to measure 
specifically the effects associated with the speed of propagation of gravity c g as it can 
be confused with the speed of light c. Thus, it is expedient to keep the notations for 
the speed of gravity and that of light different. If the retardation of gravity effect 
in the Lienard-Wiechert solution (jl.6|) of the Einstein equations could be observed it 
would allow us to determine the speed of gravity c g by testing the magnitude of the 
time-dependent terms in an affine connection and/or curvature tensor. 

1.4- Basic Principles of the Speed of Gravity Measurement by VLBI 

A simple idea for confirming the general relativistic prediction that gravity propagates 
is to measure the direction of the gravitational force exerted by moving celestial bodies 
on each other, for example, in a binary system. At the first glance this direction is 
expected to coincide with the retarded positions of the bodies on their orbits because of 
the finite speed of propagation of gravity c g ITT] . Then, experimental observation of 
the retarded positions would set up a limit on c g . Unfortunately, this method does not 
work out in a straightforward way. As a matter of fact, during the time of propagation 
of gravity between the massive bodies they move with almost constant velocities and 
acceleration-dependent terms are small. It was proved that in the case of a uniformly 
moving body the gravitational force measured at the point x is not directed from this 
point towards the retarded position of the body x a (s) (s = t—r a (s)/c s ) but points to the 
present position of the body x a (t) taken at the time of observation [TTj. Thus, measuring 
the direction of the gravitational force in the case when all acceleration-dependent terms 
are neglected does not allow to observe the effect of the propagation of gravity, that is it 
does not allow to determine whether the influence of the gravitational force comes out 
of the retarded or present position of the body. This occurs because of the aberration 
of the gravitational force in general relativity which depends on the velocity of the 
body v a and the speed of gravity c g in such a way that it compensates for the effect 
of the propagation of gravity [llj. The same kind of compensation of the retardation 
in propagation of electromagnetic field from a point charge by the aberration of the 
electric force occurs in Maxwell's electrodynamics that is a well-known phenomenon 
|12j . We shall discuss the aberration of gravity effect in relation to the jovian deflection 
experiment in section 15.31 

This cancellation between the retardation and aberration of gravity effects may lead 
to the erroneous conclusion that the gravity field does not propagate at all in the case of a 
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uniformly moving body which is equivalent to the assertion that the gravitational field 
propagates with infinite speed c g = oo [13]. Such a conclusion, however, contradicts 
the nature of the retarded Lienard-Wiechert solution of the Einstein equations and 
the causality principle of special and general relativity. Analogy with electrodynamics 
reveals that the correct statement would be that propagation of the gravitational field 
from the uniformly moving body to the field point is done in the form of the gravitational 
wave with both infinite wavelenght and period such that the speed of gravity c g remains 
constant (compare with ^3] which discusses emission of electromagnetic radiation by 
a uniformly moving charge). Purely Newtonian gravitational force with the position of 
the bodies taken at the retarded time would violate the conservation laws as first noted 
by Laplace (see, for example, for historical remarks). Additional velocity-dependent 
terms in the gravitational force must be added to compensate for this retardation in the 
Newtonian force. Thus, existence of the aberration of gravity terms in the gravitational 
force is required by the Lorentz-invariance and it implies that the gravitational field 
of the moving bodies is to be taken into account through the retarded position of the 
bodies ■ This preserves validity of the third Newton's law in relativistic binaries 

at least in conservative post-Newtonian approximations. 

In the case of a uniformly moving body the metric tensor can be represented as a 
function depending on either present x{t) or retarded position x(s) (s = t — r/c g ) of the 
body. This is because the metric tensor in the frame with respect to which the body is 
moving, can be obtained either directly in the form of the retarded Lienard-Wiechert 
solution of the Einstein equations or by means of a kinematic Lorentz transformation 
of the static (Schwarzschild) solution of the Einstein equations to the moving frame. 
If one denotes the speed of gravity c g in the Einstein equations, a matrix of the 
Lorentz transformation must be taken dependent on the speed-of-gravity parameter 
c g as well in order to transform the Schwarzschild solution from one frame to another 
without violation of the Einstein equations. Hence, the fact, that the metric tensor 
given in terms of the retarded Lienard-Wiechert potentials can be re-expressed in terms 
of the static Schwarzschild solution by making Lorentz transformation is simply the 
other form of the statement that gravity propagates with the speed c g . Had the 
matrix of the Lorentz-transformations of the Einstein equations depended on some other 
parameter c* rather than c g , would lead to mathematical inconsistency in two methods 
of obtaining of time-dependent gravitational field of a uniformly moving body. Einstein 
had equated the speed of gravity c g to c to preserve the special relativistic character 
of the Lorentz-invariance of equations of general theory of relativity. In other words, 
Einstein extrapolated his principle of special relativity to general relativity by assuming 
that the fundamental constant c characterizing the ultimate speed of propagation of 
gravitational field is the same as in special theory of relativity where gravitational field 
does not matter. This Einstein's postulate is, in fact a hypothesis which can be tested 
in gravitational experiments conducted in time-dependent gravitational fields % 

% Making experiment in the gravitational field of a uniformly moving body is sufficient. 
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We have suggested a straighforward way to test equality of c g and c [21 ITS| . 
Our method is based on measuring a gravitational perturbation of the phase of an 
electromagnetic signal coming from a quasar and passing through the time-dependent 
gravitational field of moving body (Jupiter). This perturbation of the phase is an 
integral along the light-ray trajectory which we were able to perform and express as a 
logarithmic function of a scalar product $(s) between two vectors taken at the point of 
the present position of the light particle, x a = (ct,x), where t is the current instant of 
time. A first vector k a is a wave vector of the electromagentic signal coming from the 
source of light (quasar) to an observer. A second vector r a = x a — x®(s) is directed from 
the present position of the light particle to the retarded position of the ath light-ray 
deflecting body x%(s) = (cs,x a (s)) taken on its world line at the retarded instant of 
time s defined by equation (jl.7j) that defines the retarded gravity cone along which the 
gravity field propagates from the moving body to the field point with speed c g . 

The scalar product of the two vectors is defined by equation 

$( s ) = - ka r a = r a -k-r a} (1.8) 

where k a = (1, k), k is a unit vector from the source of light to observer, r a = x — x a (s), 
and the retarded time s is given by equation (|1.7|) . Notice that in general case the 
spatial directions of the two vectors are different and do not coincide, so that ^ 0. 
It makes the light and gravity propagational directions separated in the space-time so 
that one can distinguish between them. We have found (see [2| and next sections of this 
paper) that light propagating past a massive body moving with velocity v a experiences 
a relativistic time delay 

A~^^ln$( S ), (1.9) 

where M a is mass of the light-ray deflecting body. 

The principle idea of the experiment [2] was based on the following arguments. 
One measures the relativistic time delay A and, hence, the scalar product $(s). The 
direction to the source of light k is determined by VLBI directly from measuring the 
unperturbed part of the phase of the electromagnetic wave. Direction of the vector 
r a{s) = x — x a (s) from observer to the body (Jupiter) is a function of a single parameter 
c g which magnitude fixes position of the body on its world-line with respect to the point 
of observation (t,x) via the retarded time s defined by equation (|1.7|) . By fitting the 
theoretical prediction for the scalar product, $theor y (s), to its observed value, $ bs, we 
can derive from minimization of the residual value, <5$ = | < l ) theor( s ) — "^obsl, the parameter 
c g - the (ultimate) speed of gravity - and to evaluate if it is equal to the speed of light c 
or not (see Fig. |5J). The VLBI jovian deflection experiment has proved that the ultimate 
speed of gravity is equal to the speed of light with the precision of 20% [T|. 

1.5. Notations and Conventions 

Four-dimentional coordinates on the space-time manifold are denoted as x a = (x°, x l ) = 
(ct, x l ), where c is the speed of light. Roman indices run from 1 to 3 while Greek indices 
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run from to 3. Repeated indices assume the Einstein summation rule. Indices are 
rised and lowered with full metric g a p. The Minkowski metric of flat space-time is 
Va/3 = diag(— 1, +1, +1, +1) and the Kroneker symbol 5^ = diag(l, 1, 1, 1). Partial 
derivative d a = d/dx a . Round brackets around indices denote symmetrization, for 
instance, = 1/2 (A^ + A v »). 

We also use boldface italic letters to denote spatial vectors, for instance, a = a 1 = 
(a 1 , a 2 , a 3 ). A dot between two spatial vectors is the Euclidean scalar product, i.e., 
a-b = a l b l + a 2 b 2 + a 3 b 3 . A cross between two spatial vectors is the Euclidean vector 
product, i.e., ax b = eijka^b k , where e^k is the anti-symmetric Levi-Civita symbol such 
that £123 = +1. 

2. The Speed-of-Gravity Parameterization of the Einstein Equations 

2.1. The post- Newtonian Expansion and the Speed-of-Gravity Parameterization 

Let us assume that one works in an arbitrary reference frame with coordinates x a = 
(x°,x l ). The post-Newtonian approximation scheme of solving the Einstein equations 
treats a derivative of the metric tensor g^ v with respect to time coordinate x° = ct 
as a small quantity compared with a derivative of this tensor with respect to spatial 
coordinates x % [7J |HJ E3 • I n order to parameterize this difference it is customary to put 
a small parameter e in front of each time derivative of the metric tensor explicitly, that 
is to make use of the replacement 

c dt * c dt ' 1 ' ' 

In the post-Newtonian approximation scheme the Newtonian limit of general relativity is 
achieved for e — > 0. This is equivalent to the statement that the speed of propagation of 
gravitational interaction, c g , goes to infinity (c g = oo) because in the Newtonian limit the 
gravitational field propagates instanteneously by definition [T7] . It is obvious that in this 
limit the Einstein gravity field equations must reduce (at least in one preferred frame) 
to partial differential equations of elliptic type which can not contain time derivatives 
of the gravitational field variables (the metric tensor) at all. 

From this consideration it follows that the parameterization (|2.1j) is equivalent to 
the statement that the post-Newtonian parameter e can be defined as 

e = c/c g , (2.2) 

where c is the speed of light, and c g is the speed- of- gravity parameter running from c g = c 
to c g = oo (0 < e < 1). The speed of light c pertains to all non-gravitational phenomena 
taking place in a flat space-time while the speed of gravity characterizes maximal rate 
of temporal variation of gravitational field and the speed of its propagation. For this 
reason it is natural from a physical point of view to correlate each time derivative of the 
metric tensor g a p with the speed of gravity c g . This is in accordance with approaching 
general relativity as a spin-2 field theory with self- interaction [T%1 119j. From the field- 
theoretical point of view only the speed of gravity c g is allowed to appear in the Einstein 
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tensor (see equation (|1.4jl ) because there are no other propagating fields involved in its 
construction. We shall call the post-Newtonain parameterization of the time derivatives 
of the metric tensor given by equations (j2.1|) and (|2.2j) as the speed-of-gravity (or c g ) 
parameterization of the Einstein equations. 

A four-dimentional generalization of equation (|2.1j) is achieved after introducing an 
auxilary unit vector field V a which makes the time derivative in (j2.1|) to be Lorentz- 
invariant. The field V a is adynamic and introduces preferred-frame effects which 
disappear in general relativity when c g = c. The parameterized partial derivative of 
the metric tensor g^ u in arbitrary reference frame can be written down as follows 

5 Q = d a - (e - l)V a V% = ed Q - (e - 1) P^dp , (2.3) 

where d a = d/dx a , P a ^ = i] al3 + V a V^ is the projector on the plane orthogonal to the 
four- vector V a + that allows us to preserve the Lorentz-invariance of the derivative (S a . 
In what follows, we shall simplify our equations by chosing the system of geometrized 
units such that the universal gravitational constant G and the speed of light c are equal 
to unity: G = c = 1. In this system of units the post-Newtonian speed-of-gravity 
parameter e = l/c g . 

2.2. The Einstein Field Equations 

In general relativity the metric tensor g a p relates to Ricci tensor R^ v via Christoffel 
symbols |H1 |H] 

1 

flU 



= ^ (<W + d u g p , - d p9llv ) , (2.4) 



D _ f) p" _ o -pa , -pa -pa _ -pa per /p r\ 

The speed of gravity parameterization of the Einstein equations must retain these two 
relationships but replace the derivatives d a — > B a . Thus, we obtain the parameterized 
affine connection 

f V = \ 9 aP + %9pu - 3^) , (2.6) 

and the parameterized Ricci tensor 

D _ a ~T a 3 f« I pa pc pa pwr /p y\ 

One introduces a parameterized covariant derivative for geometric sector of the 
parameterized Einstein equations in such a way that for any vector field A 11 

V a ^ = g a ^ + f^. (2.8) 



The differentiation rule (J2.8)) is compatible with the metric tensor 

V a g a ? = , (2.9) 
+ V a V a = g aP V a V^ = -1. 
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and preserves the validity of the Bianchi identity for the parameterized Einstein tensor 
Gfj, u = R^ v — (l/2)g flv R 

V a G aP = . (2.10) 

However, it violates the local equation of motion for the stress-energy tensor T a/3 in 
the sense that <5 a T al3 ^ in flat space-time. This means that the matter part of the 
Einstein equations should be generalized to make it compatible with the conservation 
laws expressed in terms of the 'crossed' post-Newtonian derivative (|2.Hj) . This goal is 
achieved with the following parameterization of the stress-energy tensor 

where the parameter 5 = e — 1. It is straightforward to prove that if d a T aP = 0, then 
3 Q G"^ = 0. Conservation law for Q al3 is now a consequence of the Bianchi identity for 
the parameterized Einstein tensor G^ v . The Einstein equations parameterized with the 
speed of gravity parameter c g reads 

GV = 8vr0^ . (2.12) 

Bianchi identity ()2.10j) assumes that 

V a Q at3 = , (2.13) 

and this equation is in concordance with the flat space-time conservation law for the 
stress-energy tensor of matter T a ^. 

We introduce the weak-field decomposition of the metric tensor 

9a/3 = Vap + Kp , (2.14) 

where h a p is the perturbation of the Minkowski metric tensor 7] a p. In what follows we 
restrict ourselves with the linear approximation of the Einstein equations with respect to 
the metric perturbations h a p. However, we shall keep all powers of the parameter e in the 
post-Newtonian expansion of the linearized theory *. In the linearized approximation 



1 

and the Ricci tensor 



R^u = \ (cTcV/w + cPcl./w - K„ - V&h) , (2.16) 
where h = rj^h^. 

Let us introduce a new field variable which is a linear combination of the metric 
tensor perturbations 

= hT - ]^rf v h . (2.17) 

The linearized gravitational field equations are obtained from the Einstein equations 
(ETrejl after making use of equations (Q, (jZHJ), (PHEJ), (^TTjl and omitting all 
quadratic and higher order perturbations in h^ u . One gets 

g°a a7 ^ _ 5^5 a7 ra - + v^^/3l al3 = -167t6^ . (2.18) 

This is called the first post-Minkowskian approximation of general relativity |15| . 
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The parameterized field equations (|2.18|) are gauge invariant with respect to a 
special group of gauge transformations of the metric perturbations that is found by a 
simple inspection of these equations. Let 7^ be an 'old' field, and 7^ is a 'new' field 
related to the 'old' one by the gauge transformation 

= y v - - g"£" + if 3^ . (2.19) 

We can check that the gauge transformation ()2.19|) does not change the field equations 
()2.18|) for the 'new' field. The gauge function £ a satisfies the homogeneous equation 

Q a 3 a £P = . (2.20) 

We shall use the gauge freedom to simplify the field equations. This will be done in the 
next section. 



2.3. The Post- Newtonian Metric Tensor 

In this section we solve the linearized field equations ([2. 18)1 . To this end let us impose 
the gauge condition 

5,7^ = , (2.21) 

on the metric tensor. This simplifies equation ()2.18j) and reduces it to 

d a Q a Y" = - 16tt0^ , (2.22) 
where the differential operator 

Q a d a = U-{e 2 -l){V a d a ) 2 

is a generalization of the Dalambert wave operator □ = r] a ^d a dp for the case of c g 7^ c. 
One can easily see that the imposed gauge condition leads to the conseravtion law for 
6^, i.e. 3^6^ = 0, that follows from equation 1|222|) . 

We concentrate on a case of an isolated N-body system like our solar system. Solar 
system moves with respect to a preferred frame with velocity V % . Preferred frame is 
often identified with isotropy of the cosmic microwave background radiation [22] • Lunar 
laser ranging proves that the preferred frame effects caused by spatial components of 
V a are small and can be neglected for analysis of the jovian deflection experiment. We 
should also keep in mind that our goal is to single out effects associated with the speed 
of gravity c g in general relativity. As soon as they are identified, we shall take a limit 
c g — > c, which will eliminate all preferred frame effects from final results irrespectively 
of the value of the spatial components of the velocity V a . 

For aforementioned reasons, we shall simplify our calculations by working in the 
reference frame of the solar system (t, x) with the origin located at its barycenter 
and assuming that in this frame the vector field V a = (1,0,0,0). It preserves all 
effects associated with c g but supresses preferred frame effects which are irrelevant for 
theoretical interpretation of the jovian deflection experiment. In the chosen frame one 
has: 5 — e<9o ; an d 9j = 9j, so that the linearized field equations ()2.22|) are simplified 

-e 2 ^ + V 2 ] Y v (t, x)=- 167r6^(t, x) , (2.23) 
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where 

@oo = T oo ^ Q oi = eT oi ^ Qlj = e 2 T io _ ( 2>24 ) 
Let T a/3 be the stress-energy tensor of massive point-like particles j2H] 

N 

T^(t,x) = Y, M «X\t)<(t)u» a {t)^\ X -x a (t)) , (2.25) 

a=l 

where the index a — 1, 2, N enumerates gravitating bodies of the solar system, M a is 
the (constant) rest mass of the ath body, x a {t) are time- dependent spatial coordinates 
of the ath body, v a (t) = dx a (t)/dt is velocity of the ath body, ti" = 3 a (l, v a ) is the 
four-velocity of the ath body, 1 a = (1 — v 2 )^ 1 ^ 2 is a Lorentz-factor, and o~( 3 )(a;) is a 
3-dimentional Dirac's delta-function. 

A physical (causal) solution of equation ()2.23|) can be found in terms of the retarded 
Lienard- Wiechert tensor potentials |20j 

7 00 (t, x) = 4 £ 1 S Ma2 H i \ > ( 2 - 26 ) 
^ r (s) - ev a (s) ■ r a (s) 

„ i( )=4f ^ A4W£M(£) , (2 . 27) 
^ r a (s) - ev a (s) ■ r a (s) 

a=l 

t^i r a (a) - ev a {s) ■ r a {s) 

that depend on the retarded time s = s(t,x) determined as a solution jj of the gravity 
cone equation 

s = t - er a (s) , (2.29) 

connecting the retarded space-time position (s,x a (s)) of ath body and the field point 
(t, x), and one uses notations r a (s) = x — x a (s), r a (s) = \r a (s)\. 

The gravity cone equation ()2.29|) describes the propagation of the gravitational field 
from the massive body to the field point [20J with the speed c g = e -1 . In the limit of 
e — > 1 the gravity cone is a null hypersurface in Minkowski space-time. In the Newtonian 
limit e — >• gravitational interaction is instanteneous (c g — > oo) and the gravity cone 
collapses to the spacelike hypersurface of constant time fT7\ . 

It is interesting to notice that parameter e = l/c g appears only in 7°* and ^ 
components of the metric tensor perturbations and it does not reveal itself in 7 00 . This 
preserves the Newtonian limit from degeneracy irrespectively of the numerical value of 
the speed of gravity c g . In other words, we keep constant the ratio of bare value of 
the universal gravitational constant G to c 2 g equal to the ratio of the observed value of 
the universal gravitational constant G to the speed of light c: G/c 2 g = G/c 2 . Thus, we 
conclude that experiments testing relativistic effects in a static field limit can not set 
up any limitation on the speed of gravity constant c g . Only gravitational experiments 
testing relativistic effects of time-dependent fields can be used to measure c g . 

t) Notice that the retarded time s is different for each body as it depends on the body's position x a . 
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One should not confuse the speed of gravity parameter c g with the parameter a\ 
of the standard PPN formalism [22] • Indeed, if one compares our equations (|2.26J1 - 
()2.28|) for the metric tensor perturbation with the corresponding components of the 
PPN metric tensor [22] one will find that the PPN parameter 7 = 1 and 



In the rest frame of the solar system and under condition 7 = 1, the parameter 
«i describes physically the coupling of the g Qi component of the metric tensor with 
matter's current, while c g defines the maximal rate of change of time derivatives of 
the gravitational potentials, that is the ultimate speed of gravity. Equation ()2.3()j) is 
a formal mathematical consequence of our intention to preserve compatibility of the 
Bianchi identity with the local laws of conservation of the stress-energy tensor of matter 
and the arrangement of the parameters in Will-Nordtvedt's PPN formalism [22 . Jovian 
deflection experiment is not sensitive enough to measure relativistic effects associated 
directly with g 0i component of the metric tensor and, hence, can not measure the 
parameter a±. But it is sensitive enough to measure contributions of the first time 
derivatives of goo and gij components of the metric tensor and confirm the Einstein 
general principle of relativity for gravitational field by measuring c g . GP-B gyroscope 
experiment [21] will be able to constrain the PPN parameter a\ in near future so that 
the consistency of equation ()2.30|) will be tested experimentally 

A standard post-Newtonian approach consists in a further expansion of equations 
(I2.2(i|) - (j2.28)l with respect to the small parameter e and dealing with each term of 
the expansion separately. However, such development of slow-motion expansion is not 
appropriate for this paper as we are looking for a calculation of the time delay in compact 
(unexpanded) form which will take into account all powers of the parameter e. This will 
allow us to compare the results of the present paper with those obtained in our preceding 
papers [21 E] and to discuss the alternative interpretations of the experiment proposed 
in [23 EH 123 EH] • An exact theory of the relativistic time delay also elucidates the 
physical origin of the retardation, observed in the orbital position of Jupiter in the VLBI 
experiment pQ, as caused by the finite time taken by gravitational field to propagate 
from moving Jupiter to the present position of a photon emitted by the quasar (see Fig. 
[0 Fig. E]). This retardation confirms validity of the Einstein principle of relativity for 
gravitational field. 

3. The Post-Newtonian Time Delay 

3.1. General Formalism 

In general relativity light propagates in vacuum along a null geodesic world line with 
the locally measurable speed c. The undisturbed propagation of light in the absence of 
gravity is then a straight line 




(2.30) 



x x N (t) = x + ck\t - t ) , 



(3.1) 
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where to is time of emission, x is coordinate of the source of light at time to, k % is 
the unit vector in the direction of propagation of light from the source to observer. A 
first approximation of the relativistic equation of light propagation parameterized by 
coordinate time t reads 

d 2 x % 



dt 2 



c 2 k^ {k» 



(3.2) 



where k^ = (1, k l ) is the four-vector of the unperturbed light-ray path. We emphasize 
that the Christoffel symbols in equation ()3.2|) depend on the parameter c g and 
are given by equation (|2.15J) . This distinguish the fundamental speed c g of the 
Einstein principle of relativity for gravitational field from the physical speed of light 
c. Christoffel symbols should be taken on the unperturbed light-ray trajectory (J3.ll) . 
that is F" = F" (t, x N (t)). Dependence of the Christoffel symbols F° on the ultimate 
speed of gravity c g is a natural property of their pure geometric origin which is not 
correlated with electromagnetic waves and the speed of their propagation. Briefly, the 
Christoffel symbols has no electromagnetic nature and can not contain characteristics 
of electromagentic field like the speed of light. Therefore, in general relativistic limit, 
when c g = c, the constant c in the Christoffel symbols must be understood as the 
ultimate speed of gravity defining in accordance to the Einstein principle of relativity 
the maximal speed of propagation of gravitational field. 

Double intergation of the ordinary differential equation ()3.2|) along the unperturbed 
light ray yields the time of propagation of light from the point x l to x\ = x l (t\) 

1, 



ti — tc 



\x 



x \ + A(t l5 t j 



where the relativistic time delay 



A(ti,t 



Wk v k c 



dt 



to 



dr 



x=x n (t) 



(3.3) 



(3.4) 



One emphasizes that in order to calculate the integral in the right side of equation (|3.4j) 
one must first take a derivative and then substitute the unperturbed light ray trajectory 
x n (t) for x ff . 

The important note here is that in accordance with our post-Newtonian 
parameterization ()2.3|) and the rule of the directional derivative 29J, one has 

d 



ck a 



where 5 



x=x N (t) 

1, and 



dr 



h 



fJtU 



J, X N {T) 



+ 5 



dh^(r,x) 
Or 



(3.5) 



- x=x n (t) 



dr dr dx % 



denotes the total time derivative along the unperturbed light ray path. By making 
use of equation (|3.5|) one recasts equation ()3.4j) for the relativistic time delay into the 

ff It is worthwhile to point out that the parameter r used in this section as a variable of integration, 
should not be confused with the parameter r used in paper [Sj. 
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following form 



A(ti,t ) 



dtl fcTL (t,x N (t)) + 5 



to 



dr k»k p 



dhjw (r, ac) 
dr 



(3.6) 



a;=ccjv(r) 



It is remarkable that the speed of gravity parameter e = c/c g appears in this equation 
not only in the metric tensor (see equations (j2.26|) - (j2.28|) ) but explicitly in front of the 
second integral. This is a crucial feature of general relativistic parameterized approach 
required for unambiguous interpretation of the relativistic VLBI experiment under 
discussion. The presence of the second integral in equation (j3.fij) was missed in papers 
of other researchers EDI- This is the mathemtical origin of the disagreement 

between our calculation of the jovian deflection experiment and that present in |2"5|l26j. 
Substitution of equations (j2.2fij) - (j2.28j) into the time delay (jH.fi |) yields 

A(ti, t ) = Ai + A 2 , 

where 



(3.7) 



dr 



1 - efe • w (s))' 



1 



a=l 

A 2 = 25^M a 

a=l 



tl 



dt 



r a (s) 



ev n (s 




x=x n (t) 



(3. 



(3.9) 



to 




1 - ek ■ v a (s)Y 



'l 



x=x n (t) 



and the gravity cone equation 

s = r — e\x — x a (s)\ , (3.10) 

gives the retarded time s = s(r, x) in the form of the implicit function of time r which is 
used in the integrands of equations (j3.8j) and (j3.9j) . The integrals in equations (j3.8j) and 
(j3.9j) have general form and are valid for arbitrary moving light-ray deflecting bodies. 
The method of calculation of these integrals have been worked out in the case of e = 1 in 
a series of our papers [HH EU EH] • The case of e ^ 1 is more involved and requires some 
simplifying assumptions. We shall suppose that the bodies deflecting light are moving 
along straight lines with constant velocities. The integrals in equations (j3.8jl and (j3.9j) 
are calculated under this assumption in the next sections. 



3.2. Primary Time Delay Integrals 

The relativistic time delay is found by performing integration in equation (j3.fi)) . When 
looking for the post-Newtonian corrections to the Shapiro time delay, positions of moving 
gravitating bodies in the metric tensor (j2.26j) - (j2.28j) are to be considered as functions 
of time x a (t). The integral (j3.fi jl can be calculated analytically for arbitrary orbital 
motion of the light-ray deflecting bodies in case of the parameter e = 1, that is c g = c 
[21 E2 EHj However, the case of e ^ 1 is too complicated and can not be treated 
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analytically, if the true (circular or elliptical) orbital motion of gravitating bodies is 
used. For this reason, we shall have to recourse to simplification of the problem and 
focus on a solvable case of each body moving with a constant speed along a straight line 
from — oo to +00. Approaching the problem in this way gives a very good approximation 
in the calculation of the relativistic time delay (and deflection of light) so far as body's 
acceleration can be neglected that is determined by the magnitude of residual terms 
and experimental accuracy. We notice that an accurate integration assumes that one 
can not replace the problem of light propagation in the field of moving bodies with 
one in which each body is fixed at the time of the closest approach of light ray to it 
|3H EE] • Such an approach leads to an approximate solution which does not preserve 
the Lorentz invariance of the theory and lead to incomplete conceptual understanding 
of the physical origin of the velocity-dependent corrections to the Shapiro time delay 

[23 ESI ESI ES!- 

Our theoretical approach assumes that trajectory of a-th body is a straight line 

x a {t) = x a {t a ) + v a (t - t a ) , (3.11) 

where v a is a constant vector of body's velocity with respect to the barycenter of the 
solar system, and t a is a fiducial instant of time that can be taken arbitrary. Time t a 
drops out of the final result and, hence, plays no role in discussion of physics of the 
experiment. This is because we consider motion of the light-ray deflecting bodies as 
uniform and rectilinear from —00 to +00. 

Let us introduce a space-like vector of a relative distance 

R a = x - x a (t) , (3.12) 

connecting the field point x a = (t, x) to the position of a-th body taken on the 
hypersurface of constant time t passing through the field point. Taking into account the 
retarded gravity cone equation ()2.29|) and expanding body's coordinate x a {t) around 
retarded instant of time, s, we can recast vector R a to another equivalent form 

R a = r a (s) - ev a r a (s) , (3.13) 

where r a (s) = x — x a (s), and r a (s) = \r a (s)\. Notice that vector r a (s) belongs to the 
surface of the gravity cone and connects the field point (t, x) with the retarded position 
of the light-ray deflecting body (see Fig. |3J). This simple transformation gives an idea 
how the aberration of gravity works to compensate the retardation of gravity effect in 
the case of a uniformly moving body. Indeed, the retardation in the orbital position 
of the body, x a (s), caused by finite speed of gravity propagation is cancelled by the 
aberration-like, velocity-dependent term in equation ()3.13j) which is also proportional 
to the speed of gravity c g . 

In the case of c g = c the vector r a (s) is a null vector, otherwise it is a space-like 
vector. Taking the square of equation (j3.13|) and solving the quadratic equation with 
respect to r a (s) yields 



ev a ■ R a ± \j R 2 a - e 2 (R a x v 
1 - e 2 v 2 



r.(«) = • ( 3 - 14 ) 
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By making use of equation ()3.13j) we can also prove that 

r a {s) - ev a ■ r a {s) = (l - e 2 v 2 a ) r a {s) - ev a ■ R a . (3.15) 

Comparing equations ()3.14j) with ()3.15j) and taking into account that the distance r a (s) 
is always positive and must coincide with R a in the limit e — > 0, we finally obtain 

r a (s) - ev a ■ r a (s) = y/l% - e 2 (R a x v a f . (3.16) 

This formula will be used for the calculation of integrals from the metric tensor (|2.26J1 - 
(12.28)) that originally depended on the combination of the retarded distances shown in 
the left side of equation ()3.16j) . 

We express the right side of equation ()3.16|) as an explicit function of time by 
substitutinging equations ()3.1|) and ()3.11|) to equation ()3.12j) . This yields the Euclidean 
distance between the light ray particle and a-th body 

R a = p(t - t ) + r , (3.17) 

where 

p = k-v a , r = x - x a (t ) , 

and the position of a-th body at the time of emission of light must be understood as 
calculated with the help of the following formula: x a (t ) = x a (t a ) +v a (t — t a ). Inserting 



equation ()3.17|) to equation ()3.15|) leads to the factorization 

R 2 a - e 2 (Ra x v a f = A-(t) ■ A+(t) , (3.18) 
where functions 

A-(t) = cq-(t - t ) + h- , (3.19) 

A+(t) =cq + (t-t ) +h+ , (3.20) 

and constant vectors in equations ()3.18j) - (|3.2()j) are defined as follows 

Q- =p-e(px v a ) , (3.21) 

q + =p + e(p x v a ) , (3.22) 

h = r - e(r x v a ) , (3.23) 

h + = ro + e(r x v a ) . (3.24) 



One can observe that the left side of equation ()3.16|) . expressed in terms of the retarded 
time s, can be written down in terms of instanteneous quantities by making use of 
equation ()3.18|) . This allows us to calculate integrals from the retarded expressions in 
the time delays given by equations (|3.7|) - (|3.9|) . 

3.3. Calculation of the Primary Integrals 

The formalism of the preceding section has been developed under the assumption that 
the light-ray deflecting bodies had constant velocities. This assumption reduces the 
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A 1 = 2G^M a Z a /«, 



0=1 

AT 



A 2 = 2G^M a Z a /( 2 ) 



a=l 
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(3.25) 
(3.26) 

(3.27) 
(3.28) 

F3/2( r ) ' (3 ' 29) 

where P = p ■ v a = k ■ v a — v\ , and Q = r ■ v a are constants, and F(t) = A_(t) ■ 
is a quadratic function of time 

F(t) = a(t-t ) 2 + /5(t-to)+7, (3-30) 
with constant coefficients 



where the Doppler factor 

= (1 -ek-v a ) 2 

a ~ ^A^! ' 

is a constant relativistic parameter. The primary integrals are 

-* 1 dt 



to 



J( 2 ) = (e-1) Trft T [P(r-t )+Q] 
</io J —00 



a = <?„ • q+ = p 2 - e 2 {p x v a ) 2 



/3 = Q ■ h+ + q+ ■ h = 2p ■ r + 2e 2 [(p • v a )(v a ■ r ) - v 2 a (p ■ r 
7 = h ■ h + = r 2 - e 2 (r x u a ) 2 . 

Performing an integration in equation ()3.28|) yields 



In 



a 



\/otF{t) 



F(t) 



t=t 



(3.31) 
(3.32) 
(3.33) 



(3.34) 



where the overdot denotes the time derivative: F(t) = dF(t)/dt, and we have shown 
explicitly the upper and lower limits of integration at which the integral must be taken. 
Double integration in equation ()3.29|) yields 

t=ti 



1 



j(2) 



(e-l)-lP ( 



1 aQ--/3P 



a 



2a 3 / 2 



F(t) 



(3.35) 



t=t 



Integrals (|3.34jl . ()3.35|1 can be expressed in terms of original variables (that is coordinates 
and velocities of the bodies and the light ray) after making use of the following exact 
relationships 



a 



1 — k ■ v n 



(fc X Vgf 

l-k- V a ) 



(3.36) 
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aQ - l -(5P = (1 - e 2 v 2 a ) [(fc • p)(« a ■ H Q ) - (k ■ R a )(p ■ «„)] , (3.37) 

= (1 - e 2 ^) (fe • «a) " (1 - e 2 fc • « a ) (v a ■ R a ) . (3.38) 

We draw the attention of the reader to the fact that the quantity aQ — (l/2)/3P in 
equation ()3.37|) is constant that can be expressed in terms of the vectors r and v a 
because all time-dependent terms in the right side of equation (J3.37)) are mutually 
canceled out. We have expressed aQ — (l/2)/3P in terms of the time-dependent distance 
R a between the light particle and a-th body as this form is more convenient for making 
further transformations. 

It is also important to notice that the product 

V^m-^j (v^ + ^j =«F(t)-^ = trr-2 ( 3 - 39 ) 

represents a constant determinant of the quadratic function F(t). This helps to 
transform original results of the integration ()3.35j) obtained in this section to a different 
form by observing that 

V^)-^ 1 = p 2 -- (3-40) 

/ a 7 - T 

This formula is not directly used in the present paper but it is important in discussion 
of the time delay formula obtained by making use of the advanced, instead of retarded, 
Lienard-Wiechert potentials. 

The primary function representing the result of the integration ( |3.34jl . (J3.35)) is 
expressed in terms of the instanteneous variables as follows 



V^W) - *~y~ = y/0- ~ k ■ v a y + ( e 2 - l)(fc x v a )*y/Rl - e?(R a x <p (3.41) 
- (1 - e 2 v 2 a ) (k ■ R a ) + (1 - e 2 k ■ v a ) (v a ■ R a ) . 

We shall use the results of this section to obtain the final expression for the 
relativistic time delay. We shall give this expression in terms of the retarded variables 
r a = x — x a (s) associated with the retarded time s given by the retarded gravity cone 
equation 1)2.29)1 . This is because such an approach fits naturally with the Minkowski 
world of special relativity and allows us to derive equations which are apparently 
invariant with respect to the Lorentz transformations in case of c g = c. Furthermore, 
the time delay in terms of the retarded variables obtained in the present paper can be 
easily compared with our previous results published in other papers. 

The Lorentz transformation technique has been used by Klioner [36J to calculate 
the relativistic time delay and deflection angle. These calculations are in agreement with 
our results which helps to understand better how the Lorentz transformations work in 
the problem under discussion. We notice also that the Lorentz transformation technique 
does not simplify calculations and/or reduce their amount. In a sense, it is technically 
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simpler to obtain the time delay and the deflection angle by direct integration of the 
light-ray equations in a moving frame. Lorentz transformation technique operates both 
with the Einstein equations and equations of light geodesies. Lorentz transformations 
bring about the time-dependent terms of the metric tensor when transforming the 
equations from static to moving frame. The time-dependent terms of the metric tensor 
are conceptually coupled with the ultimate speed of gravity. If parameter c of the 
Lorentz transformations was not equal to the ultimate speed of gravity c g one would 
observe it in the experiment under discussion. This discrepancy was not observed with 
accuracy of 20% pQ. 

3.4.. The Primary Integrals in Terms of Retarded Variables 

The primary integrals can be transformed to the retarded variables. Retarded variables 
reflect the nature of the relativistic time delay more adequate than the instantaneous 
functions. This is because our world is a four-dimentional space-time manifold with a 
causal structure determined by null hypersurfaces. This causality is apparently seen if 
one operates with the retarded variables. Moreover, equations for time delay and the 
light deflection expressed in terms of the retarded integrals can be applied to the case 
of gravitational lensing by stars of our galaxy and/or the Magellanic Clouds jHZllSHI- It 
is impossible to parameterize data processing algorithm of the microlensing events with 
the position of the lensing star taken at the present instant of time (time of observation). 
Position of the lensing star must be referred to the retarded instant of time which is 
solution of the retarded gravity cone equation (j3.10J) . 

We use equation ()3.13|) for transforming the space-like vector R a to its retarded 
countepart r a . We also remind the reader that according to its definition, y/F(t) = 
r a — et>a • r a . Thus, equation (j3.41j) can be readily recast to 



V^W) - = (r« - ev a ■ r a ) y/(l-k-v a y + (e*-l)(kxv a ) 



(3.42) 



- (l - e 2 vl) [kr a - e(k ■ v a ) r a ] + (l - e 2 k ■ v a ) (v a ■ r a - ev 2 a r a ) , 

which should be compared with equation (j3.41|) . This formula is exact and valid both 
for non-relativistic (v a < c) and relativistic (v a ~ c) velocities of the moving bodies. 
Constant term defined by equation ()3.37|) can be re-written in terms of the retarded 
quantities as follows 

<xQ~/3P = (l-e 2 v 2 a ) [^•r- a -(fc- Va )(fc-r- a ) + (fc-r a )^-(fc- Va )K-r a )- e (fcx Va ) 2 ] .(3.43) 

Now we are prepared to give the post-Newtonian expansion of the time delay and analyze 
the physical origin of various terms in this expansion. 
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4. The post-Newtonian Expansion of the Time Delay 

4-1. Time Delay and the Propagation Speed of Gravity 

The time delay is an observable quantity which is expressed in the c g -parameterized 
general relativity in the case of c g 7^ c as a sum of two integrals Ai and A2 (see equation 
()3.7|) ). The second integral A 2 in equation ()3.7j) is due to the parameterized relationship 
(j2.15j) between the Christoffel symbols and the metric tensor. This introduces the 
(ultimate) speed-of-gravity parameter e not only in the Einstein equations but in the 
equations of motion of test particles as well, hence, introducing the second integral A 2 
to the relativistic time delay (|3.7j) . This second integral makes a physically significant 
contribution to the time delay which eliminates ambiguity in the physical interpretation 
of linear v/c terms of the post-Newtonian expansion of time delay. 

By making use of equations ()3.42j) we can expand the time delay in powers of 
the post-Newtonian parameter 5 = e — 1. We shall retain in this expansion all terms 
up to first order in the velocities of the massive bodies. Terms of second and higher 
orders will be neglected because they are smaller than the present-day accuracy of VLBI 
observations. These terms of the second order could be a challenge for observations using 
future astrometric missions like GAIA and/or SIM. Therefore, their calculation is highly 
desirable and can be a matter of future work. 

The primary variable in the post-Newtonian expansion of the time delay will be a 
function 

$(s) =r a -k-r a , (4.1) 

where r a = \r a \, r a = x — x a (s) and the coordinates of the massive bodies x a (s) are 
calculated at the retarded time s = t — r a (s)/c g in accordance with the gravity cone 
equation ()2.29|) . Function $(s) is interpreted as a scalar product between two vectors in 
Minkowskii space-time. The vector k a = (l,fe) is always a null vector associated with 
propagation of light, and vector r a = (r a , r a ) is associated with propagation of gravity 
from the moving body. The latter vector is a space-like vector for any value of e 7^ 0, 
and a null vector if e = 1. 

We express the time delay in terms of function $(s). First, we have 

Fit) 



V^W) - = + 5 [i k ■ v ^a - v a ■ r a \ + O (5 v 2 a ) (4.2) 

(k ■ v a )r a -v a -r a 



1 + 6- 



+ 0(8vl) 



r a -k-r a 

The second term in square brackets is small compared with unity in any practical 
situation. Hence, we can make use of a Taylor expansion and write down the first 
primary integral in the time delay as 

/« = - -L 1„ Hs) - S (fc • Va)Ta ~ Va ' ra +0(5vl) , (4.3) 
y/a r a -k-r a 

where for the sake of simplicity we have omitted the limits of integration which can be 
easily restored every time when necessary. 
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Post-Newtonian expansion of the second primary integral entering the time delay 
yields 

if = S(k ■ v a )ip + 5 {k ' Va) l k ' r * ) ~ Va ' ra +0(6 vl) . (4.4) 

Summing up the two integrals yields a remarkable result 

J?> + If) = + *(*•"«) ln$(g) _ 5{k . Va)+0 {6vl) (45) 
V« 

Notice that the dependence on the parameter 5 of the second term in the right side of 
equation ()4.5j) is illusory: to obtain the measured effect, we must multiply equation ()4.5j) 
for body a by its mass M a , sum over all the bodies in the system in accordance with 
equations ()3.7|) . (|3.25J) . (|3.26jl . and apply the law of conservation of linear momentum 
of the gravitating system. The second term in the right side of equation (|4.5jl vanishes 
and the final formula for the time delay is 

J^GM n l-c' l (k-v a ) f ) t=tl , „ 

A( tl ,t )= / 7= rfir \Hra-k-r a )\ +0(6v 2 a ) , (4.6) 

0=1 V i_W a/ C ^ - 1 t=to 

where r a = X]y(t) — x a (s), r a = \r a \, and the retarded time 

s = t \x N (t) -x a (s)\ . (4.7) 

Notice that the general relativistic expression ()4.6j) for time delay is given in terms 
of the retarded variables and contains no linear terms proportional to the product of the 
parameter S and the body's velocity v a , as contrasted to the paper (see also [3H] ) 
where different approach has been used for the derivation of the time delay (see jS] and 
section 7 of the present paper for further discussion) . Residual terms in equation (j4.6J) 
descibe the contribution of the second-order, velocity-dependent terms to the time delay 
in the case where c g ^ c. It is clear that in the general relativistic limit, when c g — > c, 
the residual terms vanish, and the retarded time argument s in coordinates x a (s) of 
the light-ray deflecting bodies is a purely gravitational phenomenon caused by the finite 
speed of propagation of gravity c g . It has nothing to do with the physical speed of light 
c used in observations (see Fig. 



4-2. The Differential VLBI Time Delay 

The VLBI measures the time difference 

A(ti, t 2 ) = A(t 2 , t ) - A(ti, t ) , (4.8) 

where A(t 2 , to) is calculated from the equation for A(t 1; t ) after making the replacement 
ti — > t2, Si — » S2, and x\ — > x 2 . In principle, the light ray propagating from the quasar 
to the first VLBI station moves along the unit vector k\ while that propagating to the 
second VLBI station propagates along the unit vector k 2 and these two vectors, strictly 
speaking, are different. However, this difference is of order of the diurnal (geocentric) 
parallax which is neglibly small for any of the known quasars. For this reason, we assume 
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ki = k 2 = k. Detailed study of the residual terms caused by the difference between ki 
and k 2 was given in [52*] . 

Taking equation ()4.fi|) and substituting it to equation (|4.8jl yields 

" ™ ' 1 ^ 'V^j+fsr-n.Ml (4J) 



A( fl , f2 )=2^^(l + ilf.„„)l I1 

n— 1 \ g / 



r 2 a(s 2 ) + K-r 2a (s 2 ) 



a=l 

where the residual terms ~ v\ have been neglected. We have introduced in equation (|4.9j) 
the unit vector If = —A; pointing to the quasar from the barycenter of the solar system, 
and the retarded times due to the finite speed of gravity are defined in accordance with 
equation ()4.7|) as 

si = ti--|x 1 («i)-x J (s 1 )| , (4.10) 

C g 

s 2 = t 2 |x 2 (t 2 ) -Xj(s 2 )| , (4.11) 

C g 

where t\, t 2 are times of arrival of the radio signal from the quasar to the first and second 
VLBI stations respectively. The distance between each telescope and Jupiter are given 
by r u = \r u \, r 2J = \r 2 j\, r 1J (s 1 ) = xi(ti) - Xj(si) and r 2J (s 2 ) = x 2 (t 2 ) - x,/(s 2 ) , all 
of which depend implicitly on the value of the ultimate speed of gravity c g . 

The relativistic VLBI experiment on September 8, 2002 was designed to set an 
upper limit on the speed of gravity c g by measuring the differential time delay of radio 
wave propagating through the gravitational field of moving Jupiter although the time 
delays A ffi and A caused by the gravitational fields of Earth and Sun respectively, 
are also important and had been taken into account. Expanding retarded coordinate 
of Jupiter xj(s) in equation (J4.9)) in a Taylor series around the time of observation t\ 
(that is equivalent to the post-Newtonian expansion with respect to the small parameter 
Vj/c g , where Vj is the orbital velocity of Jupiter) yields [UI2] 

A(t u t 2 ) = A ffi + A + Aj , (4.12) 

Aj = -olj-B , (4.13) 
AGM 



1 _ 2N- Vj \ N + Vj - (K-vj)K 



c g © y e 



+ 0(|), (4.14) 



c 3 Ru 

where the subscript J refers to Jupiter, Ru = x\ — Xj(ti), Ru = \Ru\, © = 
arccos(— -Ru ■ K / Ru) — |£|/ r ©J is the (small) angle between the undisturbed 
astrometric position of the quasar and the present position Xj(ti) of Jupiter at the 
time of observation, £ is the vector of the light-ray impact parameter with respect to 
the present position of Jupiter, and B = x 2 {t\) — X\{t\) is a baseline between two VLBI 
stations. The unit vector N = £/|£| is orthogonal to the unit vector K and is defined 
in such a way that the unit vector L = Ru/Ru is decomposed in accordance with 
equation [2] 

L = -if cos + AT sin9 . (4.15) 
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This equation makes it clear that the impact parameter £ = lies in the plane of two 
vectors K and L, is orthogonal to vector K (that is lies in the plane of the sky), and 
is directed from the present position of Jupiter Xj(t\) towards the light-ray trajectory. 

Formula (|4.14j) defines the angle of relativistic deflection of light as a function of 
time ti and demonstrates that there are two terms in the post-Newtonian expansion 
of the time delay caused by Jupiter. The first term describes relativistic deflection of 
light in the plane of the sky which is proportional to the static (Einstein) deflection and 
directed radially outward of Jupiter. The second term describes relativistic deflection 
of light in the direction of the barycentric velocity of Jupiter Vj taken at time t% and 
projected onto the plane of the sky. The experiment, as it was designed, was most 
sensitive to the second term because the part of the first (radial) term proportional to 
c g vanishes at the time of the closest approach t# of Jupiter to the quasar (in the plane 
of the sky) because at this time the scalar product N(t#)-Vj = 0. 

We used formula ()4.9|) in the analysis of observational data and measuring the 
ultimate speed of gravity c g |T. Formula (j4.14j) was used for making numerical 
estimates of the retardation of gravity effect and for better physical understanding of 
the relativistic deflection of light by moving Jupiter. Some authors [2Ej have stated 
that the post-Newtonian expansion ()4.14j) is illegitimate. This might be true if the 
light-ray deflecting body was very far away so that the post-Newtonian series could 
not be convergent (this is the case of microlensing events jHHEB])- However, this post- 
Newtonian expansion is rapidly convergent in the solar system and, thus, mathematically 
is well-defined because the parameter of the expansion, that is the ratio of vj/c s to the 
angle G, is limited by the value < taken on the date of the minimal approach of 
Jupiter to the quasar jlj. This value is rapidly decreasing as the angle 6 is growing. 
Hence, all conditions being necessary to apply the Taylor expansion are satisfied and 
the post-Newtonian equation (j4.14j) is fully legitimate. 

5. General Relativistic Interpretation of the Jovian Deflection Experiment 

5.1. Confirmation of the Einstein Principle of Relativity for Gravitational Field and 
Measurement of the Ultimate Speed of Gravity 

We have constructed a speed-of-gravity parameterization of the field equations of 
general relativity to separate relativistic effects associated with fundamental constants 
pertaining to gravity and electromagnetic field. We calculated a higher-order velocity- 
dependent contributions to the Shapiro time delay shown in section El Its post- 
Newtonian expansion yields formula ()4.6|) leading to the differential VLBI time delay 
given by equation ()4.9|) in terms of the retarded position of Jupiter due to the finite 
time gravity takes in order to propagate from moving Jupiter to the point of observation 
(see Fig. [Q. Expansion of the retarded arguments around the time of arrival of the 
electromagentic signal to the first VLBI station brings about equations (|4.12jl - (|4.14jl 
which are identical to the time delay equation derived in our previous works PQI2]- The 
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present paper proves that the measured parameter 5 is defined as 5 = c/c g — 1 and, hence, 
describes a difference between the speed of light c and that of gravity c g 1J. Parameter 
5 has a real physical meaning measuring discrepancy between the fundamental speed of 
the Einstein principle of relativity for gravitational field (the ultimate speed of gravity 
c g ) and the speed of light c. Measurement of this parameter in the relativistic VLBI 
experiment conducted on September 8, 2002 allows us to measure the ultimate speed 
of gravity c g . We have determined c g = (1.06 ± 0.21)c [I] and confirmed validity of 
Einstein's general principle of relativity for gravitational field. 

Deflection of light by the time- dependent gravitational fields has been analyzed in 
|HU 1321 ES] in case of c g = c. In the present paper we have established that the post- 
Newtonian formula for the time delay parameterized by the speed of gravity c g is given 
by equation 1)3. 6 j) . It is more instructive to recast it to the following form 

A(*i,*o) = ~ / dt f k»k 

^ Jt J -oo 

This equation has two terms - the first describes contribution of the space derivatives 
of the metric tensor to the time delay and the second one illustrates the contribution 
of the partial time derivative of the metric tensor to the time delay. Space derivatives 
measure a spatial inhomogeniety of the gravitational field of Jupiter as light moves 
across it. Time derivatives of the metric tensor multiplied with the parameter e = c/c g 
measure how fast a temporal change of the gravitational field of Jupiter produced by 
its orbital motion is transmitted from Jupiter to the light-ray trajectory. The rate of 
this transmission is given in terms of the dimensionless coefficient 1/e = c g /c. In the 
Newtonian gravity the speed of gravity c g = oo, the rate of the transmission is infinite, 
and there is no contribution of the time derivative of the metric tensor to the integrated 
time delay of light. Physically it means that the Einstein principle of relativity for 
gravitational field is violated and gravity propagates so fast that the gravitational field 
has no time to change during the time of flight of light ray. On the other hand, in the 
case of c g = c the partial time derivative of the metric tensor gives a contribution to 
the time delay but it is smaller than that of the spatial derivatives by the factor v a /c 
as expected in general relativity which assumes that the Einstein principle of relativity 
is valid for gravitational field. Despite the smallness of the time-dependent effect our 
experiment has convincingly demonstrated that it exists and definitely contributes to the 
integrated time delay because the speed of gravity c g = (1.06 ± 0.21)c is not infinite. In 
other words, the jovian deflection experiment confirms the Einstein principle of relativity 
for gravitational field with accuracy 20%. 

5.2. Testing the Weak Principle of Equivalence in Time- Dependent Gravitational Field 

The weak principle of equivalence states that test bodies move along geodesies of curved 
space-time which locally coincide with straight lines of Minkowski geometry j?!] . So far, 
the weak principle of equivalence was tested only in static gravitational field of Earth 
and Sun [22] but it was not known whether it is valid in time- dependent gravitational 



dx l 



+ e 



dh^ (t,x) 
dr 



dr 



(5.1) 
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fields or not. The quantitative measure of validity of the weak principle of equivalence in 
time-dependent gravitational field is parameter e = c/c g which characterizes magnitude 
of the contribution of the time derivatives of metric tensor to the affine connection. 
According to general relativity parameter e = 1 which means that the metric tensor is 
fully compatible with the affine connection in the sense that the covariant derivative of 
the metric tensor is zero irrespectively of whether gravitational field is static or time- 
dependent. Jovian deflection experiment measures c g and consequently provides a test 
of the weak principle of equivalence for photons in time-dependent gravitational field. 
Our measurement of c g proves that with the precision of 20% the weak principle of 
equivalence for time-dependent gravitational field is valid and the fundamental constant 
c g coupled with the time derivatives of the metric tensor in the affine connection 
coefficients is numerically equal to the speed of light c. 



5.3. Measuring the Aberration of Gravity Field 

Let us introduce two angles G = ©(£) and 9 = 9(s) between the vector k characterizing 
direction of propagation of the light ray and the unit vectors L = R\j / Ru and 
/ = rij/ru, where Ru — X\ — xj{ti) connects present position of Jupiter Xj(ti) 
at the time of observation t\ and the point of observation X\, and ru = x\ — Xj{s\) 
connects the retarded position of Jupiter at the retarded time s\ = t\— ru/c g and the 
point x\. By definition cos = k ■ L and cos 9 = k ■ I (see Fig. QJ. One recalls that we 
have already used definition of the angle O in equation (|4.14)1 . 

The logarithm in the time delay equation ()4.6|) is ln(ru — k ■ ru) = In ru + ln(l — 
cos#), so that the general relativistic time delay A j(ti, to) due to Jupiter can be written 
as 



Aj(ti,t ) 



2GM J l-c g 1 (k- Vj ) 



In ru + ln(l — cos i 



(5.2) 



c 3 y/l - v'j/c* 

where we have dropped the constant term at the time t of emission of light by the 
quasar. The retarded coordinate of Jupiter xj(si) can be expanded around the time 
t\ in a Taylor (post-Newtonian) series with respect to the time difference s — t with 
subsequent replacement of this difference by making use of equation (|2.29|) . Accounting 
for that linis-^ 9 = O it gives 

9 = e-(r 1 j/c g )e + 0(r 2 u /c 2 g ), (5.3) 

where 

• k ■ v j — (k ■ L)(L ■ v A 

Ru sin 6 

The difference between the two angles <fi = 9 — O is called the aberration of gravity angle. 
This is because it is this angle which determines the difference in the direction of the 
Newtonian gravitational force of Jupiter calculated with or without accounting for finite 
speed of propagation of gravity and measured at the point of observation by making 
use of photons propagating in the time-dependent gravitational field. The aberration of 
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gravity effect is elusive and can not be observed for slowly moving bodies at the order 
of (f/c) 3 [TH I13j. However, it can be observed for photons because they move with 
the same speed as the speed of gravity propagation. The aberration of gravity effect 
would disappear if the speed of gravity c g = oo because in such a case 9 = O. From the 
definition of the angle 9 and equation ()5.3j) one can represent the aberration of gravity 
effect in vector form 



ln(l - cos#) = ln(l -kl)= In 



k ■ ( L + —L x (vj x L 

Ca 



0[ V 1). (5.5) 



It is important to realize that in the VLBI experiment under discussion the only 
object directly observed in the sky was the quasar, that is we measured components of 
the vector k very precisely. Jupiter was not observed directly in radio. This is because 
there is virtually no emission that is detected from Jupiter with VLBI observations. 
Any radiating object more than about 10 mas in size would be invisible. Even if a 
few percent of the emission from Jupiter did come through the VLBI observations, one 
would have no idea of where is was coming from and it probably would be a combination 
of small features on Jupiter which happen to have enhanced emission. We implicitly 
measured the position of the center of mass of Jupiter by our detection of the Shapiro 
delay. Our accuracy of the Shapiro delay was about lknow where Jupiter was, if we 
knew its mass we could deduce its position to an accuracy of 3.7/100 or about 2" 
(two arcseconds) accuracy. However, the JPL ephemerides for Jupiter are much more 
accurate than this, and this is the position that we used P] . For this reason, neither the 
angle 9 nor the angle were directly measured. The value of was precisly determined 
at the time of observation from JPL ephemerides while the retarded angle 9 in the sky 
with respect to the quasar was obtained by fitting theoretical model of the relativistic 
deflection of light from the quasar to observations. This part of the experiment was 
drastically misunderstood by Samuel |2E : who assumed that we measured position of 
quasar with respect to Jupiter by measuring relative position of the quasar with respect 
to Jupiter in radio. This led him to the wrong conclusion that effects of order vj/c 
beyond the Shapiro delay could not be observed (see section [7| for discussion of this 
Samuel's misinterpretation in more detail).. 

Our experiment measures the magnitude of the aberration of gravity effect which 
is inversly proportional to the speed of gravity c g from Einstein's general principle of 
relativity. Measuring the aberration of gravity allows to determine the speed of gravity c g 
in the jovian deflection experiment in close analogy with the Bradley's idea of measuring 
the speed of light from observation of the aberration of light. The aberration of gravity 
effect can be also re-formulated in terms of gravitodynamic dragging of the light ray from 
the quasar caused by the translational motion of gravitational field of Jupiter jB] . It can 
be viewed as a generalization of the Fizeau effect for the case of moving gravitational 
field considered as a medium in flat space-time with an effective refractive index defined 
by the Newtonian gravitational potential. 
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6. Relativistic Time Delay in Two-Parametric Model of Gravity 

While the speed-of-gravity parameterization of the Einstein equations of general 
relativity can be done with a single parameter c g , one can admit existence of alternative 
theories of gravity which may have different values of the ultimate speed of gravity 
parameter, c g , and the speed of gravitational waves parameter, cgw. In this section we 
shall assume that general relativistic relationship between the Christoffel symbols and 
Ricci tensor is violated and the vacuum equations of an alternative gravity field theory 
depend on two parameters - c gw and c g such that c gw 7^ c g (see section IT~T|) . There 
exists no yet a convincing example of a self-consistent theory of gravity of this type. 
Hence, in order to proceed we shall simply corrupt relationships of Einstein's general 
relativity to get some insight to those effects which can be observed in the light-deflection 
type experiments in the field of moving bodies. Such corruption however can not be 
extrapolated too far because any corruption of fundamental laws of gravitational physics 
reduces our ability for adequate physical interpretation of observed facts. For example, 
a simple geometric picture of the effect measured in the jovian deflection experiment 
and shown in Fig. (JHJ does not applicable anymore in the two-parametric model of 
the gravity "theory" and we are lost in the midst of zillion possible interpretations 
while only one is supposed to be true in nature. To restrict the freedom existing in the 
two-parametric model we shall preserve the Einstein principle of relativity according to 
which the speed of gravitational waves cgw can not exceed the ultimate speed of gravity 
c g that limits the rate of change of time derivatives of the metric tensor. Indeed, it is 
highly unlikely from the point of view of modern physics to expect that gravitational 
waves would propagate with arbitrary speed exceeding the ultimate speed of gravity c g 
which in its own turn is to be equal numerically to the speed of light. 

In the two-parameteric model of gravitational equations the parameter c g is now 
enters into all time derivatives of the metric tensor excluding the second time derivative 
entering the wave operator and depending by definition on the other constant c gw . We 
continue to denote the partial time derivatives with c g as 5 a = (3o,9i) = (c" 1 ^,^). 
The partial derivatives with c gw will be denoted as d a = (d , c\) = (c~^d t , di). 

In the linearized approximation we shall have by definition: 

• the affine connection 

• the linearized "Einstein" tensor 

qv» = _ d*d a Y» + 5^7^ + d"d a Y a ~ r/^SaS/37 Q/3 , (6.2) 
where 7^ is defined by equation (j2.17|) ; 

• the linearized "Einstein" equations 

d a a ^ U - ^Q a Y a - S"Sa7^ + = - 16tt0^ . (6.3) 

The field equations ()6.3|) are not gauge-invariant and, already for this reason, can not 
represent a self-satisfactory model of an alternative theory of gravity. In principle, 
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additional terms must be added to the left side of equation (|6.3|) to recover the gauge 
invariance of the theory. It requires much more work and it is very likely that the 
result will be negative, that is the gauge-invariant field equations of the two-parametric 
"theory" do not exist. For this reason we use equations (|6.3|) just formally without 
giving them too much physical content. This agrees with the point of view expressed 
in [3U] who used another alternative two-parametric model of the field equations and 
showed that interpretation of the relativistic time delay in such model crucially depends 
on the relationship between the speed of light c and the speed c gw . 

The law of conservation of the stress-energy tensor demands that the following 
differential conditions on the field variables have to be imposed 

5V7^ = . (6.4) 

Then, the field equations ()6.3|) are reduced to the following form 

d a d al ^ = - 16tt0^ . (6.5) 

Condition ()6.4|) restricts the dynamical freedom in chosing gravitational variables. 
However, one must keep in mind that it is imposed by the law of conservation of matter 
and does not reflect a real dynamic gauge freedom of the gravitational field variables 
which does not exist in this model. 

Repeating the calculations of the time delay as presented in section 01 one obtains 
(compare with equation (|4.6|) ) 
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One can observe that in the two-parametric model of the gravity field 'theory' 
(which is too unrestrictive because of the non-gauge invariance) two speed parameters 
can be independently measured in time delay - Cgw and c g . However, a scrutiny 
examination reveals that the parameter shows up only in the second order terms of 
the post-Newtonian expansion of equation (|6.6j) . Linear terms of the post-Newtonian 
expansion of equation ()6.6)1 contain only the ultimate speed of gravity c g . Therefore, 
without having extra information about the relationship between c g and the 
measurement of the linear v/c terms beyond the Shapiro delay can not tell us anything 
about the speed c gw without making use of an additional information. This information 



CONTENTS 



34 



is provided by the Einstein general principle of relativity for gravitational field which 
says that as soon as the ultimate speed of gravity c g is known the speed of gravitational 
waves Cgw < c g for any viable alternative theory of gravity. Equality c gw = c g is achieved 
for massless gravitons in general relativity. Hence, the jovian deflection experiment sets 
an upper limit on the speed of propagation of gravitational waves c gw irrespectively of 
the model one uses for fitting the observations. Current precision of the jovian deflection 
experiment sets the gravity speed limit c g = (1.06 ± 0.21)c [T]. 

One-parametric general relativistic approach used in the present paper predicts 
that = c g and this is what we should expect when future space missions like 
GAIA dhttp://astro.estec.esa.nl/GAIA/] ), SIM flhttp: / /planetquest.jpl.nasa.gov/SIM/| ) , 
or ARISE [^n] will presumably allow the measurement of quadratic corrections to 
the Shapiro time delay and/or light deflection angle and, thus, disentangle the speed 
c gw from c g directly. Such relativistic experiments will provide us with much deeper 
knowledge on the structure of time-dependent terms in equations of the gravitational 
theories and, for this reason, they are highly desirable for inclusion to the working plans 
of these and future space missions. 

We would like to emphasize that the long-term observations of the binary pulsar 
PSR 1913+16 [S] alone do not allow us to establish that c g = Cgw- This is because if one 
introduces two parameters c g and Cgw separately to the timing model then the number 
of unknown fitting parameters will exceed the number of the observed post-Keplerian 
parameters. Hence, the system of equations for finding all fitting parameters can not be 
resolved. However, if one makes use of the limit on c g set up by the VLBI experiment 
under discussion jT], an upper limit on Cgw from the binary pulsar data is Cgw < 1.27c. 
Relativistic binary pulsars with large proper motion provides acess to larger number of 
observable parameters jlUj and they may be used for separate measurement of c g and 
Cgw without resorting to any other kind of observations. 

7. Alternative Interpretations of the Jovian Deflection Experiment 

Difficulties in the interpretation of the jovian deflection experiment reflects a hidden 
variety of physical meanings of the fundamental constant c entering differing physical 
equations. Usually this constant is called in relativistic sleng as "the speed of light" 
without paying respect to the nature of the equations it enters in. This may be 
misleading in the case of the affine connection and curvature tensor entering equations 
of general relativity since it can provoke a wrong understanding of the real nature of 
these objects which is not electromagnetic but purely geometric. Affine connection 
and curvature characterize the intrinsic properties of the space-time manifold which are 
described by general relativity being conceptually independent of the Maxwell theory. 
Dynamical properties of the geometric sector of general relativity are characterized by 
the fundamental constant c but it has physical meaning of the ultimate speed of gravity 
[4~T] . General relativity provides a tool for separation of the propagation of gravity 
effects from those associated with light. This tool is the Lienard-Wiechert solution of 
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the Einstein equations which connects the field point with the gravity-generating body 
by a null cone establishing causal character of the propagation of the gravitational field 
due to the finite speed of gravity. Solution of Maxwell's equations in the gravitational 
field does not affect the structure of this gravity null cone which is a primary object. 
Phase of electromagnetic field is perturbed by the moving gravitating body from its 
retarded position defined by the solution of the retarded time equation which is an 
integral part of the gravitational Lienard-Wiechert potentials. Hence, the retardation 
in position of the light-ray deflecting body is a unique property of gravity pointing out 
to its causal nature and finite speed of propagation. The fundamental constant c from 
the geometric sector of general relativity (affine connection, curvature tensor) can not be 
measured in laboratory. Only gravitational experiments in time-dependent gravitational 
field can provide us with its numerical value. Jovian deflection experiment is the first one 
where the fundamental speed c from the geometric sector of general relativity has been 
measured and the Einstein general principle of relativity was confirmed for gravitational 
field. 

Some researchers [211 E3 EH1 12H] claimed that the jovian deflection experiment is 
sensitive only to the physical speed of light that is the fundamental constant c, which 
we measured has an electromagnetic origin associated with Maxwell equations. In other 
words, these researchers disagree that the jovian deflection experiment measures the 
fundamental speed of the Einstein general principle of relativity for gravitational field 
(the ultimate speed c g of gravitational interaction). Let us discuss their arguments. 

The interpretation of the experiment by Asada (2Z| is that the retarded position 
of Jupiter measured in the experiment is due to the propagation of the radio waves 
from the quasar. He noticed that the observer, source of light (quasar), and Jupiter 
are lying on a surface of a past null cone with vertex at the point of observation (see 
FigEJ). The light ray of the quasar moves to observer along the hypersurface of this 
null cone. According to Asada [21], Jupiter can be considered as fixed in space at the 
retarded position Xj(s) during the time of the experiment. A null line connecting the 
retarded position of Jupiter and the observer belongs to the hypersurface of the past 
null cone. Because a radio wave from the quasar moves along the light cone and the null 
line connecting Jupiter and the observer also belongs to this past null cone, Asada has 
concluded that this must be the physical speed of the radio wave c entering equation 
(II. 7|) instead of the speed of gravity c g . Thus, the nature of the retardation in position 
of Jupiter was explained in [21] as due to the finite speed of light from the quasar. 
However, physical light does not propagate along the null line connecting Jupiter and 
observer. This null direction is accociated with the causal character of the gravity field 
equations and its observation in the experiment tells us that the speed of gravity is 
the same as the physical speed of light with accuracy 20% 1 . A correct Minkowski 
diagram of the experiment is shown in Fig. and includes two null cones one of which 
describes causality region of the gravity field from Jupiter and the second one - does 
causality region of light from quasar. Additional discussion of Asada's interpretation 
can be found in j^]. 
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Samuel [2E] made use of somewhat convoluted considerations partially involving 
the Lorentz-transformation technique and derived the time delay in terms of the 
retarded angle 9 = 9(s) defined in the present paper in section 15. HI and shown in 
Fig. HJ The final formula of Samuel's paper |2E] for the differential time delay 
£2) = A(i 2 , to) — A(ii, t ) is given by his equation (8). We use our precise equation 
(|5.2|) to compare with equation (8) of Samuel's paper and to separate effects associated 
with the speed of gravity from those correlated with light. To this end we shall asssume 
that the angle 9 is small so that one can use approximations sin# = 9 + 0(9 3 ), and 
cos 6* = 1 — (l/2)9 2 + 0(9 4 ). We shall also make use of the fact that the VLBI baseline 
B = X2(t±) — Xi(ti) between two VLBI stations is small compared with the distance 
ru = \xi — xj(s)\. It allows us to expand function A(t 2 ,t ) in Taylor series with 
respect to the baseline B and the VLBI time lag t 2 — t\ = (K-B)/c, where c is the 
speed of light. One should take into account that the retarded time s 2 = t 2 — r 2 j(s 2 )/c g , 
where r 2J = \r 2J \, r 2J = x 2 (t 2 ) — Xj(s 2 ), is not a constant but a function of the point 
x 2 (t 2 ) = x 2 (ti) + v 2 (ti)(t 2 — ti) + with v 2 being the barycentric velocity of the second 
VLBI station. Differentiation of the retarded time equation ()2.29j) yields 

?L=-L r A . (7.1) 

dx l c g r j - c~ l vj ■ rj 

If one neglects the geocentric rotation of Earth, then v 2 = ve, where ve is the orbital 
motion of Earth and the Taylor expansion of r 2 j(s 2 ) provides 

r2j{s 2 ) =r 1 j{s 1 ) + B + ^^-v 2 + ^-v J , (7.2) 

c c g 

r 2 j(s 2 ) = r u (s 2 ) + I B + -(K-B)(l-v E ) + -(l-B)(bvj) , (7.3) 

c c g 

where ru(si) = Xi(ti) — xj(si), vector I = ru /\r\j\ is defined in section 15.31 and 
directed from the retarded position of Jupiter to Earth, and vector K = — k is directed 
from the solar system barycenter to quasar. 

One can decompose the unit vector I in the plane of two other unit vectors K and 
n as follows 

I = -Kcos9 + nsin9 , (7.4) 

where the angle 9 = 9(s) was defined in section 15.31 the unit vector n = n(s) = 
£( s )/l£(s)l and £(s) is the retarded impact parameter of light ray with respect to the 
retarded position of Jupiter. Vector decomposition ()7.4j) should be compared with 
equation ()4. 15|) defined at the time of observation t±. We emphasize that signs in 
the decomposition of vector I correspond to the convention of signs adopted in equation 
(j4.15|) and also correspond to that used by Samuel [2S] . The product r\j9 ~ |£(s)| = £(s) 
is equal to the magnitude of the impact parameter of the light ray with respect to the 
retarded position of Jupiter Xj(s) as defined in Samuel's paper [28j. 
Further straightforward but tedious calculations yield 
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Formula (|7.5|) coincides with equation (8) of Samuel's paper if c g = c, and we notice that 
Samuel's angle 9 ^ s coincides with our angle 9 if c g = c. The reader must also take into 
account that velocity of Earth ve and that of Jupiter vj are taken with respect to the 
barycenter of the solar system while in Samuel's paper the velocity of Jupiter (denoted 
by the letter vj) is taken with respect to Earth's geocenter. It corresponds to our time 
delay equation ()7.5|) if c g = c. Thus, we conclude that the time delay equation (8) 
from Samuel's paper is correct in the case of c g = c. However, our c g -parameterization 
of the Einstein equations reveals that Samuel has misinterpreted the meaning of the 
retarded time argument s = t — rj(s)/c g of the angle 9(s) in his derivation because he 
worked under assumption that c g = c which makes it difficult to track down the effects 
associated with the speed of gravity c g . 

Equation ()7.5)1 does not contain terms being quadratic in 1/9. It may make an 
impression that the orbital motion of Jupiter does not provide any significant deviation 
from the Einstein's prediction of the light deflection because all velocity-dependent 
terms in the right side of equation (|7.5jl are smaller than the main term (proportional to 
(n-B)/9) by a factor of 10~ 4 and can not be observed with the present-day technology 
This was the reason for Samuel's statement j2H] that terms of order v/c beyond the 
Shapiro time delay are not observable. This statement is erroneous because the Shapiro 
time delay must be calculated in terms of the present position of Jupiter at the time 
of observation t\. It corresponds to presentation of equation (|7.5jl in terms of the 
unit vector N and the angle defined in sections 14.21 and 15.31 These quantities were 
calculated on the basis of JPL ephemerides with the precision which is more than enough 
to discriminate between the time delay given in terms of the retarded variables n(s) and 
6(s) and that given in terms of the instantaneous N and G. In other words, the Shapiro 
time delay is defined as a first term in the post-Newtonian expansion of the main term 
in the right side of equation ([7.5)1 . The goal of the jovian deflection experiment was to 
distinguish two angles 9 and G. Confirmation that the apparent position of the quasar 
in the sky makes the angle 9 rather than G with respect to Jupiter is a proof that gravity 
propagates with the speed c g . 

We shall show that equation ()7.5|) is reduced to equation (j4.13[) after making its 
post-Newtonian expansion, that is the expansion of all quantities depending on the 
retarded time s\ around the time of observation t\. It reveals the linear terms of order 
v/c beyond the Shapiro time delay which can be (and have been) measured in our 
experiment [1\. We have 

-—\vj-(l-vj)l\ . (7.6) 
dsi ru L J 

Differentiation of equation ([7.4)) with respect to the retarded time Si and taking the 

small angle approximation yields 

u^^lH-ne). (7.7) 



ds\ 9 

Post-Newtonian expansion of the unit vector n = n(si) is n = N + h(si — t\) + .... By 



CONTENTS 



38 



making use of the retarded time equation s\ = t\ — r%j/c g and equation (|7.7j) one has 
n = N -^(i-nd) . (7.8) 

Futhermore, the post-Newtonian expansion of the angle 9 = 9(si) = + 9(si — t\) + ... 
after making use of the retarded time equation, reads 

e = Q-9 r -^- . (7.9) 

c g 

Combining equations ()7.8)1 and ()7.9j) we get the post-Newtonian expansion of the ratio 
n , 20 7^ AT nj Z 




\ 9 c B I c g 2 

Now we take into account that in the first approximation r±j ~ Ru, — 0, and 
9 ~ where is defined by equation ([5.4)1 . We should also use the approximation 
L ~ —K + iV0 which follows directly from equation ()4.15|) in the small angle 
approximation. This finally yields 

Plugging this formula in equation (|7.5j) brings it to the same form as the differential 
VLBI time delay given by equations (|4.13|) . ()4.14|) which was derived in our previous 
papers [21 G] arid confirmed by the experiment with the precision of 20% pQ. 

Interpretation of the jovian deflection experiment given by Will [22] (see also 26 J 
is based on the two-parametric model explained in section El of the present paper. Will 
assumed that the ultimate speed of gravity c g must be equal to the speed of light in 
accordance with the weak principle of equivalence. Hence, he called the constant c g as 
the speed of light. However, this constant is not associated with electromagnetic field 
and it is misleading to call it as the speed of light without precise explanation of its 
nature. Theoretical postulate of c g = c is that is accepted in general relativity on the 
basis of the postulate that Einstein's general principle of relativity is valid and gravity 
propagates with the same speed as light but this postulate was never tested before the 
jovian deflection experiment. 
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Figure 1. Jupiter moves in the plane of the picture (plane of the sky) from 
right to left. Positions of Jupiter taken at the time of observations t%, (2,..., £5 are 
shown by dashed circumferences and numbered as 1,2,..., 5. Retarded positions of 
Jupiter calculated at retarded times Si — ti — rj(si)/c g (i=l,2,...,5) are numbered as 
1*,2*, ...,5*. If gravity propagates with c g ~ 00 the apparent position of the quasar 
in the sky moves counterclockwise through points 1,2,. ..,5 on the green circle. If 
gravity propagates with the speed c g = c the apparent position of the quasar in the 
sky moves counterclockwise through points 1*,2*,...,5* on the red circle. The red 
circle is obtained from the green one by dragging each point of the green circle along 
Jupiter's direction of motion at the distance ~ (v.j/c g )Q~ 2 shown in equation (|7.11() . 
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Jupiter's velocity vector 




Figure 2. Absolute value of the impact parameter of the light ray with respect to 
Jupiter on September 8, 2002 is shown along with the overall magnitude of the light 
deflection (diameter of green circle) and the aberration of gravity effect (red oval) . 
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Figure 3. Light and gravity null cones are shown. They are directed to the future and 
describe propagation of light from the quasar and that of the gravity field of Jupiter 
respectively. The equation of the light null cone depends on the speed of light c while 
the equation for the gravity null cone depends on the speed of gravity c g . In general 
relativity c g = c. Both cones intersects at the field point x a = (t,x). Dashed arrow 
points from the field point to the position of Jupiter on the hypersurface of constant 
time t which was excluded by the experiment. 
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Figure 4. Light propagates from the quasar towards VLBI station (observer) on the 
Earth. While the light propagates Jupiter is moving. General relativity predicts 012] 
that the light observed at the time t, is deflected the most strongly by Jupiter when it 
is located at the retarded position xj(s) (s = t — r/c s ) irrespectively of the direction 
of propagation of the light ray and the magnitude of the light-ray impact parameter 
with respect to Jupiter. 
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Figure 5. Two null cones related to the experiment are shown. The gravity cone is a 
retarded Lineard-Wiechert solution of the Einstein equations and describes propagation 
of gravity. The light cone shows propagation of light from the quasar. The relativistic 
perturbation of a light ray measured by an observer takes place when the gravity cone of 
Jupiter passes through the observer. The VLBI experiment measures the Minkowski 
dot product $ between two null vectors at the point of observation directed to the 
quasar and to Jupiter respectively. Had Jupiter not been detected at the retarded 
position on its world line the speed of gravity c g were not equal to the speed of light 
c, and the Einstein principle of relativity for gravitational field would be violated. 
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Observer 



'Asada's "Light" Cone 




uasar 



Jupiter' s World Line 



Figure 6. Asada's light cone intersects with the quasar and Jupiter at its retarded 
position. However, no physical signals are propagating along the surface of this "light" 
cone except for two null lines which are parts of two future null cones associated with 
propagation of gravity from Jupiter and light from the quasar (compare with Fig. [5J. 
The null line connecting Jupiter and observer is a part of the gravity null cone and 
physical light from the quasar does not propagate along this line. 



